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摘 要 : 当前 ， 能 够 实现 作物 表 型 参数 高 效 、 准 确 的 测量 和 作物 生育 期 表 型 参数 的 动态 量化 研究 是 表 型 研 


究 和 育种 中 或 待 解决 的 问题 之 一 。 本 研究 以 棉花 为 研究 对 象 ， 采 用 三 维 激 光 扫 描 LiDAR 技 术 获 取 槐 人 花 植株 
的 多 时 序 点 云 数 据 ， 针 对 棉花 植株 主干 的 几何 特性 ， 利 用 随机 抽样 一 致 算法 (RANSAC) 结合 直线 模型 完成 
主干 提取 ， 并 对 剩余 的 点 云 进行 区 域 增 长 聚 类 ， 实 现 各 叶片 的 分 割 ; 在 此 基础 上 ， 完 成 植株 体积 、 株 高 、 


叶 长 、 叶 宽 等 性 状 参数 的 舍 计 。 针 对 多 时 序 棉花 激光 点 云 数 据 ， 采 用 匈牙利 算法 完成 相 邻 时 序 作 物 点 云 数 
据 的 对 齐 、 叶 片 器 官 对 应 关系 的 建立 。 同 时 ， 对 各 植株 表 型 参数 动态 变化 过 程 进行 了 量化 。 本 研究 针对 3 株 
棉花 的 4 个 生长 点 的 点 云 数 据 ， 分 别 完成 了 主干 提取 、 叶 片 分 割 ， 以 及 表 型 参数 测量 和 动态 量化 。 试 验 结果 


表明 ， 本 研究 所 采用 的 主干 提取 及 叶片 分 割 方法 能 够 实现 棉花 的 校 干 和 叶片 分 割 。 提 取 的 株 高 、 叶 长 "P 
宽 等 表 型 参数 与 人 工 测量 值 的 决定 系数 均 趋 近 于 1.0; 同时 ， 本 研究 实现 了 棉花 表 型 参数 的 动态 量化 过 程 ， 


为 三 维 表 型 技术 的 实现 提供 了 一 种 有 效 的 方法 。 
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1 引言 


近年 来 ， 随 着 农作物 表 型 组 学 和 表 型 技术 的 
发 展 ， 三 维 表 型 技术 由 于 能 够 获取 比 二 维 表 型 技 
术 多 一 个 维度 的 信息 ， 逐渐 引起 了 高 通 量 表 型 研 
究 人 员 的 关注 和 兴趣 。 其 中 ， 激 光 雷 达 (Light 
Detection and Ranging, LiDAR) 技术 作为 一 种 
通过 发 射 激光 脉冲 快速 获取 目标 空间 数据 的 有 效 
手段 ， 既 能 无 损 地 获取 目标 对 象 的 三 维 信息 ， 也 
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克服 了 传统 的 计算 机 视觉 技术 重建 作物 的 算法 复 
杂 度 高 和 耗 时 长 等 不 利 因素 ， 为 三 维 空间 中 表 型 
参数 的 自动 测量 提供 了 一 种 新 思路 。 农 作物 是 一 
种 典型 的 形态 结构 随 着 时 间 变 化 而 不 断 发 生变 化 
的 对 象 ， 对 不 同 生 育 期 的 作物 进行 持续 的 点 云 数 
据 采 集 ， 既 能 够 获取 作物 某 一 生长 点 的 表 型 参 
数 ， 也 可 以 实现 作物 表 型 参数 的 动态 变化 过 程 ， 
为 育种 学 家 提供 更 具 参 考 价 值 的 表 型 数据 。 

现 阶段 ， 作 物 的 三 维 点 云 信息 主要 通过 2 种 
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途径 获取 。 一 是 利用 基于 图 像 的 计算 机 视觉 方法 
HE (EWAt, Wwa Bin? AMOK 
米 VAR LARA) 等 。 基 于 图 像 的 计算 机 视 
觉 方 法 一 般 通 过 相机 标定 、 特 征 提取 、 特 征 匹 配 
以 及 光束 法 平 差 等 关键 步骤 实现 ， 对 获取 的 图 像 
有 约束 性 要 求 ， 如 保证 相 邻 视点 影像 间 较 高 的 重 
苹 度 ， 且 生成 的 点 云 数 据 质量 和 图 像 数 量 关 系 较 
大 。 当 作物 本 身 缺 乏 足 够 特征 时 ， 需 引入 结构 光 
等 获得 较 丰 富 的 纹理 特征 ， 以 期 获得 较 高 精度 的 
点 云 数 据 。Nguyen 等 “搭建 结构 光 系 统 获 取 植 
物 多 视 影像 ， 获 得 多 视点 云 ， 并 完成 点 云 间 的 配 
准 ， 构 建 了 植株 的 完整 三 维 模型 ; 在 此 基础 上 ， 
完成 叶片 数量 、 株 高 、 叶 片 大 小 和 节 间 距 等 参数 
的 测量 。Ni 等 ”利用 三 台 相 机 获取 120 幅 蓝莓 复 
图 像 数 据 ， 重 建 蓝 每 点 云 及 模型 。 针 对 获取 的 二 
维 图 像 数据 ， 利 用 Mask-RCNN 实 现 单 粒 蓝 侮 的 
分 割 及 成 熟 度 检 测 ， 并 将 分 割 结果 反 投 回 三 维 空 
间 ， 在 三 维 空间 中 完成 各 蓝 葡 表 型 参数 的 提取 。 
二 是 利用 LiDAR 技术 直接 获取 目标 三 维 信息 。 
该 类 技术 首先 获取 作物 对 象 不 同 视点 的 三 维 点 云 
数据 ， 再 利用 最 临近 点 迭代 法 (Iterative Closest 
Point，ICP) 完成 多 视点 的 点 云 配 准 ， 生 成 目标 
完整 的 三 维 模型 。 利 用 LiDAR 技术 开展 三 维 表 
型 研究 可 概括 为 2 个 方面 : 点 云 数据 处 理 方法 和 
表 型 参数 提取 方法 研究 。 从 作物 尺度 进行 划分 ， 
LIDAR 技术 可 分 别 应 用 在 单 株 作物 和 大 田 作物 
中 ， 对 应 地 获取 不 同 斥 度 的 表 型 参数 。 针 对 大 田 
作物 ， 提 取 的 表 型 参数 主要 包括 平均 冠 层 高 度 、 
投影 面积 、 植 株 密度 、 生 物 量 等 信息 Jiang 
E'S) 对 田间 获取 的 激光 点 云 数据 研发 了 一 套 和 
蓝 等 收获 相关 的 表 型 参数 获 的 三 维 表 型 技术 ， 包 
括 大 小 、 尺 度 等 。Jin 等 提出 了 一 种 均值 归 一 
化 向 量 生 长 算法 (Median Normalized Vector 
Growth, MNVG) 实现 田间 玉米 的 茎 叶 分 割 ， 完 
成 了 叶 面 角 、 蔡 高 、 株 高 等 表 型 参数 的 测量 。 
Sun £F "" 针对 重建 的 田间 成 熟 期 棉花 植株 ， 利 
用 密度 聚 类 方法 实现 了 棉铃 球 的 分 制 和 计数 ， 为 
棉花 估 产 提供 数据 支持 。 刘 守 阳 等 … 通过 构建 


的 数字 化 植物 表 型 平台 (Digital Plant Phenotyp- 
ing Platform, D3P) 获取 田间 小 麦 的 三 维 冠 层 场 
景 ， 记 录 了 绿色 叶 面 积 指数 、 平 均 倾 角 和 散射 光 
截获 率 等 信息 ， 在 此 基础 上 模拟 虚拟 LiDAR 实 
验 ， 生 成 了 三 维 点 云 数据 ， 展 示 了 LiDAR 数据 
在 田间 冠 层 光 截 获 和 和 冠 层 结构 方面 的 潜力 。 针 对 
单 株 作物 ， 主 要 采用 近 距 离 或 手持 式 激光 扫描 对 
象 ， 从 叶 、 茎 、 其 至 根 等 器 官 角 度 完成 各 表 型 参 
数 的 提取 。 其 中 ， 点 云 配 准 和 点 云 分 割 分 别 是 获 
取 三 维 点 云 模 型 和 提取 表 型 参数 的 关键 步骤 。 马 
福 峰 等 "提出 了 一 种 基于 植株 点 云 邻 域 几何 特 
征 约束 的 配 准 方法 ,选取 关键 点 并 估计 其 中 每 个 
点 的 支撑 邻 域 估算 出 邻 域 几何 特征 ， 利 用 几何 特 
征 实现 点 云 的 初始 配 准 。Chaudhury 和 Barron ^ 
将 不 同 视点 的 点 云 投影 至 二 维 平 面 上 ， 通 过 定位 
节点 位 置 来 建立 点 云 间 的 初始 对 应 关系 ， 该 算法 
实现 的 前 提 是 点 云 数据 中 存在 直线 特征 。 点 云 分 
割 是 获取 表 型 参数 的 前 提 。 适 用 于 作物 器 官 的 点 
云 分 割 方法 主要 3 种 : 中 基于 三 维 区 域 增长 的 方 
法 ， 通 过 利用 区 域内 部 点 云 的 相似 属性 ， 制 定 合 
适 的 区 域 增长 规则 来 完成 ， 主 要 应 用 在 植株 叶片 
EA) op EH, aA ae fF r 
等 ; @ 基 于 模型 拟 合 的 点 云 分 割 方法 ， 利 用 植株 
中 的 直线 、 圆 柱 体 、 平 面 等 典型 几何 模型 来 完 
成 ， 在 作物 主 茎 的 分 割 中 应 用 较 多 ， 如 癌 日 
3k ARK EA 等 ;加 基于 机 器 学 习 的 方 
法 ， 代 表 性 的 算法 包括 文 持 向 量 机 (Support 
Vector Machine, SVM) 和 均值 漂移 (Mean 
Shift) 等 各 种 聚 类 算法 。Wahabzada 等 ”在 获 
取 植 株 模型 的 基础 上 ， 使 用 无 监督 聚 类 的 方法 ， 
完成 了 葡萄 、 小 麦 和 大 麦 的 叶片 和 葵 秆 的 分 割 。 
喻 查 慎 等 "” 融合 了 点 云 数据 的 多 维特 征 ， 利 用 
人 工 标记 数据 为 样本 ， 结 合 SVM 分 类 器 实现 了 
树木 枝叶 分 离 ， 属 于 有 监督 分 类 方法 。 

按照 一 定 的 时 间 间 隔 多 次 采集 作物 对 象 的 点 
云 数 据 即 可 得 到 作物 的 多 时 序 点 云 数据 ， 这 些 数 
据 中 包含 了 作物 对 象 的 时 空 形态 变化 信息 。 针 对 
多 时 序 作 物 点 云 数据 的 研究 可 为 作物 表 型 参数 的 
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动态 变化 提供 强 有 力 的 手段 和 方法 ， 同 时， 也 可 
为 构建 作物 的 生长 模型 提供 数据 基础 。Su 等 2 
使 用 激光 扫描 仪 获取 了 玉米 6 个 生长 阶段 的 点 云 
数据 ， 计 算 了 玉米 株 高 、 叶 面积 和 投影 叶 面 积 3 
个 表 型 参数 ， 并 人 研究 表 型 参数 的 动态 变化 趋势 。 
Chaudhury ^j: °°") 构建 了 基于 激光 扫描 仪 的 近景 视 
觉 系 统 ， 对 野生 拟 南 芥 和 大 麦 植株 完成 不 同 视 点 
的 点 云 完成 配 准 ， 通 过 构建 表面 不 规则 网 实现 植 
株 基 于 网 格 的 面积 计算 和 体积 计算 ， 并 绘制 了 植 
株 从 观测 第 1 天 到 第 22 天 的 面积 和 体积 连续 变化 
曲线 。Sun 等 ” 使 用 激光 扫描 系统 对 大 田 棉 花 
完成 了 全 生育 期 的 表 型 参数 提取 和 动态 变化 分 
Wr, 探讨 了 表 型 参数 和 产量 之 间 的 关系 。An 
等 ” 对 拟 南 芥 进 行 了 连续 10 天 的 二 维 影像 数据 
采集 ， 分 别 从 二 维 影像 和 三 维 模型 的 角度 对 比 叶 
片 表 型 参数 测量 结果 的 差异 性 ， 结 果 表 明 与 传统 
的 二 维 影像 相 比 ， 三 维 模型 能 够 获得 精度 更 高 的 
表 型 参数 。 以 上 研究 表明 ， 虽 然 三 维 表 型 技术 已 
经 取得 了 一 定 的 进展 ， 但 针对 不 同 作物 植株 ， 数 
据 处 理 方法 各 异 ， 因 此 ， 其 潜力 还 有 待 进一步 
挖掘 。 

棉花 是 中 国 四 大 经 济 作物 之 一 ， 是 纺织 工业 
的 主要 原材料 ， 最 大 程度 地 提高 棉花 产量 是 栽培 
和 育种 的 首要 目标 ， 对 中 国 的 农业 经 济 发 展 意义 
重大 ， 而 各 生育 期 的 表 型 参数 可 为 该 目标 的 实现 
提供 数据 基础 。 为 进一步 挖掘 三 维 表 型 技术 在 单 
株 棉花 植株 上 的 淤 力 ， 本 研究 以 单 株 棉花 为 研究 
对 象 ， 利 用 激光 扫描 技术 获取 多 时 序 点 云 数据 ， 
研究 主干 和 叶片 的 分 割 算法 ， 实 现 株 高 、 叶 长 、 
叶 宽 、 叶 面积 ， 以 及 棉花 植株 体积 等 表 型 参数 的 
计算 。 在 此 基础 上 ， 实 现 棉花 植株 表 型 参数 的 动 
态 量化 过 程 ， 从 而 为 三 维 空间 的 高 通 量 表 型 组 学 
及 构建 作物 生长 模型 提供 方法 支撑 和 数据 依据 。 


2 材料 与 方法 


2.1 点 云 数据 获取 
本 研究 所 使 用 ZGScan 手 持 式 激光 扫描 仪 


(武汉 中 观 自动 化 科技 有 限 公 司 ) 获取 3 株 棉花 
共 12 个 时 间 点 的 点 云 数据 。 在 扫描 过 程 中 使 用 
的 辅助 扫描 工具 包括 激光 扫描 标记 点 、 支 架 以 及 
扫描 数据 输入 和 存储 的 惠普 笔记 本 电脑 。 使 用 戴 
尔 计算 机 进行 后 期 点 云 数 据 分 析 和 处 理 。 该 计算 
配置 为 : Windows 64 位 操作 系统 ，16 GB 运行 内 
存 ， 处 理 器 为 主 频 3.4 GHz 的 Inter (R) Core i7- 
3770 CPU。 图 1 分 别 为 移植 后 的 棉花 植株 和 扫描 
设备 。 棉 花 植株 于 2018 年 8 月 10 日 在 实验 室 中 
进行 培养 ， 待 长 至 约 16 cm 后 移 至 栽培 盆 中 
种 植 。 


(a) 移 栽 后 的 棉花 植株 (b)3D 激 光 扫 描 设备 
图 1 移 栽 后 的 棉花 植株 和 3D 激 光 扫 描 设备 


Fig. 1 The transplanted cotton plant and 3D laser scanner 

株 高 是 作物 形态 学 重要 参数 之 一 ， 叶 片 是 作 
物 进行 光合 作用 和 呼吸 作用 的 主要 载体 。 棉 花生 
长 前 期 ， 叶 和 葵 秆 的 生长 最 活跃 。 葵 秆 的 生长 ， 
叶片 数 和 叶 面 积 的 增长 为 后 期 的 营养 生长 和 生殖 
生长 积累 大 量 的 有 机 物 ， 是 影响 棉花 高 产 的 重要 
因素 之 一 。 鉴 于 葵 、 叶 对 作物 生长 发 育 的 重要 
性 ， 本 研究 主要 探讨 LiDAR 技术 在 棉花 株 高 、 
叶 长 、 叶 面积 ， 以 及 体积 等 重要 表 型 参数 测量 方 
面 的 潜力 ， 针 对 3 株 棉花 的 4 个 时 间 点 进行 点 云 
数据 采集 ， 具 体 为 种 植 之 后 的 第 46 天 、 第 49 天 、 
第 58 天 和 第 63 天 ， 共 采集 12 组 棉花 点 云 数 据 。 
随 着 棉花 植株 的 生长 ， 其 形态 结构 发 生 了 变化 ， 
主要 表现 为 植株 渐 高 ， 部 分 新 叶 萌 出 ， 部 分 叶片 
面积 增 大 ， 以 及 少数 老 叶 脱落 ， 叶 倾角 渐 大 ， 节 
间距 有 少许 增 大 等 趋势 。 使 用 Geomagic Studio 
2013 软件 对 采集 的 原始 棉花 点 云 数 据 进行 预 处 
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理 ， 包 括 删除 栽培 盆 、 文 架 等 噪声 数据 。 完 成 棉 
花 植 株 点 云 数据 预 处 理 后 ， 棉 花 植 株 的 点 云 数据 
分 别 如 图 2 所 示 。 


MRT 


(a) 棉 花 植 株 1 


ae 


(b) 棉 花 植 株 2 


(c) 棉 花 植 株 3 
注 : 每 横 排 从 左 至 右 分 别 为 第 46 天 49 天 、58 天 和 第 63 天 植株 
图 2 3 株 棉 花 的 4 个 生长 点 的 点 云 数据 
Fig. 2 Point clouds data of three cotton plants at four 


different growth stages 


2.2 棉花 器 官 分 割 方法 


后 续 人 研究 需要 通过 作物 点 云 数 据 提取 包括 株 
高 、 叶 长 、 叶 宽 、 以 及 体积 等 一 系列 表 型 参数 ， 
并 对 其 动态 变化 进行 量化 。 因 此 ， 需 要 对 棉花 植 
株 点 云 数 据 进行 器 官 分 制 ， 即 将 完整 的 棉花 植株 
点 云 数据 分 割 成 知 干 互 不 相交 的 点 云 子 集 。 
2.2.1 主干 提取 

由 图 2 可 以 看 出 ,棉花 植株 的 主干 趋 于 笔 
直 。 因 此 ,根据 棉花 植株 的 形态 结构 以 及 点 云 数 
据 特 点 ,采用 基于 模型 的 检测 方法 一 一 随机 采样 
一 致 性 (Random Sample Consensus, RANSAC) 
算法 结合 直线 模型 提取 主干 部 分 。RANSAC 算 
法 是 一 种 通过 随机 采样 观测 数据 估计 模型 参数 的 


方法 。 首 先 ， 在 原始 数据 集中 随机 指定 一 组 数据 
集合 ， 设 定 一 个 适合 该 数据 的 参数 模型 以 及 置信 
度 参 数 。 然 后 ， 通 过 迭代 选择 原始 数据 的 随机 子 
集 来 不 断 优 化 模型 的 参数 。 达 到 人 迭代 终止 条 件 
后 ， 选 取 局 内 点 数量 最 多 的 模型 即 为 所 求 结 果 。 
对 于 多 分 枝 株 型 的 作物 ， 可 以 采取 多 次 使 用 
RANSAC 的 方式 ， 直 到 所 有 的 类 直线 型 枝 干 都 
被 提取 出 为 止 。 
2.2.2 ”叶片 分 割 

在 使 用 RANSAC 方 法 对 植株 进行 主干 分 割 
后 ， 得 到 包括 多 个 叶片 的 彼此 分 离 的 点 云 数 据 ， 
只 需要 将 叶片 点 云 进行 聚 类 即 可 。 因 此 ， 对 提取 
主干 后 的 剩余 点 云 数据 使 用 欧式 聚 类 分 割 算法 即 
可 得 到 独立 的 叶片 。 欧 式 聚 类 分 割 算 法 的 步骤 
如 下 。 

(1) 为 输入 的 点 云 数据 集 P 创 建 K-D 树 ; 

(2) 构建 一 个 聚 类 C 的 空 列表 以 及 需要 检验 
的 点 的 队列 Q; 

(3) 在 当前 点 云 数 据 集 P 中 随机 选取 点 忆 ， 
执行 下 面 的 步骤。 

中 将 忆 添 加 到 当前 队列 Q; 

@ 对 于 每 一 个 属于 Q 中 的 点 己 : 

在 半径 r< qd 的 范围 内 搜索 与 点 P, 相 邻 的 点 
的 集合 Pi; 

(3 对 于 每 一 个 属于 PP 集合 中 的 点 ,检查 当 
前 点 是 否 已 经 被 处 理 过 ， 如 果 没 有 则 将 其 添加 到 
BAZ QE; 

@ 当 队列 Q 中 的 所 有 点 全 部 处 理 后 ,将 队列 
Q 中 的 点 添加 到 聚 类 C 中 ， 同 时 将 队列 Q 重 置 
为 空 ; 

(4) 当 P 中 的 所 有 点 都 被 处 理 完 并 且 属 于 聚 
类 C 中 的 一 部 分 时 结束 算法 。 
2.3 相 邻 生育 期 作物 点 云 配 准 及 叶片 对 应 

关系 建立 

为 研究 同一 株 作物 生长 变化 及 相 邻 时 期 表 型 
参数 动态 变化 ， 需 要 将 相 邻 时 期 的 作物 点 云 数 据 
DEAT ACME, ERATE PSR: (1) 点 云 数 
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据 粗 配 准 ; (2) 建立 作物 叶片 器 官 间 的 一 一 对 应 
关系 。 
2.3.1 相 邻 生育 期 点 云 粗 配 准 

点 云 配 准 根据 精度 可 分 为 粗 配 准 和 精 配 准 。 
粗 配 准 是 源 点 云 与 目标 点 云 在 完全 未 知 任何 初始 
相对 位 置 情况 下 的 配 准 。 该 操作 的 主要 目的 是 在 
初始 条 件 未 知 的 情况 下 ， 快 速 估 算 一 个 大 致 的 点 
云 配 准 和 矩阵 。 较 典型 的 是 基于 全 局 搜索 策略 的 代 
表 算 法 采样 一 致 性 初始 配 准 (SAmple Consensus 
Initial Alignment，SAC-IA) 算法 。 该 算法 在 源 
点 云 与 目标 点 云 之 间 ， 随 机 选取 几何 特征 一 致 的 
点 组 成 点 对 而 非 所 有 存在 对 应 关系 的 组 合 ， 通 过 
计算 对 应 点 对 的 变换 关系 得 到 最 优 解 ”，。 

首先 ， 计 算 点 云 的 快速 点 特征 直方 图 (Fast 
Point Feature Histograms，FPFH)。 通 过 参数 化 查 
询 点 与 邻 域 点 之 间 的 空间 差异 ， 形 成 描述 查询 点 
的 K 邻 域 几何 属性 的 多 维 直 方 图 。 

计算 每 个 点 的 FPFH 后 ， 使 用 采样 一 致 性 初 
台 配 准 算法 完成 点 云 粗 配 准 ， 具 体 如 下 。 

(1) 从 待 配 准点 云 S 中 选取 n 个 采样 点 ,为 
保证 所 采样 的 点 具有 不 同 的 FPFH 特征 ， 采 样 点 
两 两 之 间 的 距离 应 满足 大 于 预先 给 定 最 小 距离 阔 
{Ei d; 

(2) 在 目标 点 云 了 中 查找 与 点 云 S 中 采样 点 
具有 相似 FPFH 特征 的 一 个 或 多 个 点 ， 从 这 些 相 
似 点 中 随机 选取 一 个 点 作为 点 云 $ 在 目标 点 云 及 
中 的 一 一 对 应 点 ; 

(3) 计算 对 应 点 之 间 的 变换 矩阵 ， 通 过 求解 
对 应 点 变换 后 的 “距离 误差 和 ” 李 数 来 判断 当前 
配 准 变换 的 性 能 。 此 处 的 距离 误差 和 水 数 多 使 用 
Huber 误 差 水 数 表示 ， 记 为 五 (1,))， 见 公式 (1)。 


In 四 «m, 
H(1)- : (1) 
z(l4l- m). Ii >m, 
其 中 ， 肥 为 预先 给 定 值 ; /为 第 ;组 对 应 点 变 
换 之 后 的 距离 差 。 上 述 配 准 的 最 终 目的 是 在 所 有 
变换 中 找到 一 组 最 优 解 ， 使 得 误差 函数 的 值 最 
小 ， 此 时 的 变换 即 为 最 终 的 配 准 变换 矩阵 ， 进 一 


步 可 得 到 配 准 结果 。 
2.3.2 ” 相 邻 生育 期 棉花 叶片 对 应 

为 量化 作物 表 型 参数 的 动态 变化 过 程 ， 在 作 
物 点 云 数 据 粗 对 齐 的 基础 上 ， 需 要 建立 植株 叶片 
间 的 对 应 关系 。 

对 于 同一 株 作物 也 和 也 ,时 期 的 两 组 点 云 ， 
作物 7 时 期 叶片 数 为 m，7,,, 时 期 叶片 数 为 n。 
随 着 植物 的 生长 ， 可 能 会 发 生 新 叶 长 出 或 老 叶 脱 
落 等 现象 ， 因 此 ， 叶 数 m 和 nn 可 能 不 同 ,下面 将 
分 别 讨论 叶片 数 m 和 nn 相同 和 不 同 的 两 种 数据 处 
理 情况 。 


注 : 从 棉花 点 云 距 离 阅 值 左 至 右 分 别 是 1.0、5.0、6.0 和 18.0 mm 
(a) 棉 花 点 云 在 不 同 距离 阔 值 下 的 主干 提取 结果 


^. 7 


79.75. 
o~ ac fa 
| | | | 


注 : 第 从 左 到 右 分 别 为 待 分 割 叶片 点 云 最 小 聚 类 点 数 和 距离 阅 值 


min min 


nan=800 , t5120 mm 


(b) 待 分 割 叶片 点 云 和 不 同 最 小 聚 类 点 数 和 距离 闪 值 的 分 割 效果 


图 3 棉花 主干 提取 和 叶片 聚 类 结果 


Fig. 3 Results ofmain stem extraction and leaf clustering 


(1) 当 m=n 时 。 计 算 7, 和 7,, ,时 期 植株 每 个 
叶片 的 中 心 点 ， 然 后 分 别 计算 T 时 期 每 个 叶片 中 
心 点 到 7,, ,时 期 每 个 叶片 中 心 点 的 三 维 欧 氏 距 
BP, OKTAT, ,时 期 叶片 中 心 点 的 三 维 欧式 距离 
构建 一 个 m x n 的 邻接 和 矩阵 MM，M Gj) 表示 也 
时 期 第 i 个 叶片 中 心 到 工 , ,时 期 第 个 叶片 中 心 的 
三 维 欧 式 距离 。 设 有 2 组 叶片 分 别 为 L,、L,、 工 ， 
和 ZL'、L,、L;， 图 3 描述 了 这 两 组 叶片 各 中 心 点 
的 欧式 距离 。 和 矩阵 的 行 标 和 列 标 分 别 代表 7 和 
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T AE A io TA T, ,时 期 植株 的 点 
去 数据 处 于 同一 个 坐标 系 下 ， 使 用 匈牙利 算法 
(Hungarian Algorithm) XJ AE PE M 求 解 ， 使 得 
相 邻 时 期 叶片 中 心 点 的 三 维 欧式 距离 之 和 最 小 ， 
此 时 即 为 相 邻 时 期 植株 叶片 间 的 最 佳 匹 配 ， 如 
表 1 中 加 粗 距 离 值 所 示 。 可 得 出 结论 ，L, 和 工 ,为 
同一 组 叶片 ，L, 和 工 ' 为 同一 组 叶片 ， 而 工 则 和 工 ; 
形成 配对 关系 。 
表 1 不 同时 间 点 同一 植株 叶片 欧式 距离 矩阵 构建 


Table 1 Distance matrix among all the leaves of one cotton 


plant at two different time points 


TAREE Fr 
了 时 叶片 
Ly Lj L, 
Li 3.0 2.0 4.0 
L, 1.0 2.0 3.0 
"m 2.5 3.0 2.0 


(2) 4Sm<nitt. XF E — PRTEVI T, I T, 
两 个 时 期 植株 有 效 叶 片 数目 不 同 的 情况 ， 将 也 时 
期 的 叶片 数目 补充 为 n 个 ， 同 时 将 7 时 期 n 一 m 
个 叶片 中 心 到 7, ,时 期 每 个 叶片 中 心 的 距离 都 记 
为 0， 并 构建 一 个 n x n 的 矩阵， 然后 使 用 匈牙利 
算法 来 对 该 矩 阵 进 行 求解 ， 使 得 7, 和 7,, ,两 时 期 
叶片 中 心 三 维 欧式 距离 之 和 最 小 ， 此 时 得 到 7 时 
期 区 个 叶片 与 7,, 时 期 mw 个 叶片 的 对 应 关系 。 


也 ,时 期 剩余 2 一 天 个 叶片 在 到 时 期 无 叶片 与 其 
对 应 ， 则 记 为 新 长 出 的 叶片 。 

各 生育 期 的 棉花 叶片 间 一 一 对 应 关系 建立 
后 ， 即 可 完成 株 高 、 叶 长 、 叶 宽 等 一 系列 表 型 参 
数 的 计算 ， 并 可 以 构建 各 表 型 参数 的 时 序 变 
化 图 。 


3 试验 和 结果 分 析 


3.1 棉花 主干 和 叶片 分 割 结果 


RANSAC 实 现 作 物 主干 提取 的 理想 的 结 
是 叶柄 与 主干 分 离 ， 主 干部 分 的 噪声 少 并 且 主 干 
以 外 的 部 分 叶片 数 与 人 工 计 数 结果 相同 。 在 使 用 
RANSAC 算 法 实现 作物 主干 提取 过 程 中 ,不 同 
的 参数 对 结果 产生 不 同 的 影响 。 点 到 模型 的 距离 
国 值 过 小 ， 则 提取 到 的 主干 不 人 全， 距离 浆 值 过 大 
则 不 能 正确 地 提取 主干 部 分 。 如 果 和 迭代 次 数 过 
少 ， 则 不 能 得 到 最 优 的 模型 局 内 点 ， 迭 代 次 数 过 
大 会 增加 算法 运行 时 间 。 因 此 ， 合 适 的 参数 选择 
至 关 重 要 。 在 迭代 次 数 固定 的 情况 下 ， 分 别 讨论 
T EX 1.0. 6.0 Fl 18.0 mm 的 主干 提取 结 
如 图 3 (a) 所 示 。 经 过 多 次 试验 发 现 当 迭代 次 数 
为 2000 次 ,距离 阔 值 设置 在 5.0—10.0 mm 时 ， 
如 表 2 所 示 ， 可 以 得 到 较为 理想 的 棉花 植株 主干 
提取 结果 。 


表 2 棉花 株 高 测量 结果 
Table 2 The experimental results of plant height 


结果 评定 1 号 植株 


2 号 植株 3 号 植株 


时 间 / 天 46 49 58 63 46 
人 工 量 测 /mm 364.60 371.10 460.80 520.10 391.50 
株 高 计算 值 /mm 362.50 369.50 459.30 517.90 389.20 
绝对 误差 /mm 2.10 1.60 1.50 2.20 2.30 
相对 误差 /% 0.58 0.43 0.33 0.42 0.59 


49 58 63 46 49 58 63 
426.50 515.60 578.20 312.80 344.80 415.60 480.60 
424.30 513.40 574.80 31040 341.60 413.30 477.90 

2.20 2.20 3.40 2.40 3.20 2.30 2.70 

0.52 0.43 0.59 0.77 0.93 0.55 0.56 


对 于 多 分 枝 株 型 的 作物 ， 可 以 采取 多 次 使 用 
RANSAC 的 方式 ， 直 到 所 有 的 类 直线 型 枝 干 都 
被 提取 出 为 止 。 使 用 RANSAC 算 法 进行 分 割 的 
评价 标准 根据 目的 而 定 。 在 本 试验 中 ， 目 视 评 价 
指标 是 主干 应 能 反映 出 作物 的 株 高 ; 量化 
RANSAC 算法 提取 主干 精度 的 评价 指标 是 主干 


实际 长 度 与 提取 主干 长 度 的 比值 ， 比 值 越 接近 1, 
精度 相对 越 高 。 

欧式 聚 类 分 割 算 法 的 理想 结果 是 分 制 后 的 聚 
类 的 叶片 数量 与 人 工 计数 结果 相同 且 分 割 结果 中 
叶片 完整 。 算 法 实现 棉花 叶片 的 分 割 过 程 中 ， 需 
要 设置 三 个 参数 ， 近 邻 搜索 半径 、 最 小 聚 类 点 数 
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和 最 大 聚 类 点 数 。 最 小 聚 类 点 数 和 最 大 聚 类 点 数 
用 于 排除 噪声 以 及 避免 输出 近邻 搜索 半径 设置 不 
佳 时 的 错误 结果 ， 根 据 叶子 的 点 云 分 布 情况 ， 一 
般 最 大 聚 类 点 设置 大 于 点 云 数 量 的 一 半 、 最 小 聚 
类 点 设置 不 超过 点 云 数量 的 10%; 近邻 搜索 半径 
和 点 与 点 之 间 的 距离 均值 有 关 ， 如 果 近 邻 搜索 半 
径 过 大 ， 则 点 云 无 法 分 制 ， 如 果 过 小 则 会 导致 分 
制 结果 错误 。 经 过 一 系列 的 试验 验证 ， 各 参数 如 
下 设置 时 ， 可 以 取得 较 好 的 分 类 结果 : 最 小 聚 类 
点 数 为 800 个 ， 最 大 聚 类 点 数 为 135,000 个 ， 近 
邻 搜索 半径 为 120 mm。 如 图 3 (b) 所 示 ， 最 小 
聚 类 点 数 为 100 个 时 ， 主 干 提取 后 留 下 的 非 叶片 
点 云 噪声 被 单独 分 割 出 来 ， 对 最 终结 果 产 生 影 
响 ; FEES REY 80 mm 时 ， 由 于 叶片 区 域 平 均 点 
间距 大 于 叶柄 区 域 ， 阔 值 小 于 叶片 区 域 平均 点 云 
间距 时 ， 叶 片 被 分 割 为 大 量 细 碎 的 点 ， 因 此 在 最 
小 聚 类 点 数 的 限制 下 只 有 较为 密集 的 叶柄 部 分 和 
叶片 中 的 密集 区 域 会 被 输出 。 故 当 提 高 最 小 聚 类 
点 数 到 800 个 后 ， 非 叶片 点 云 噪 声 被 排除 ; 而 提 
高 距离 靖 值 到 120 mm 后 ， 叶 片 与 叶柄 便 可 以 构 
成 一 个 整体 。 

在 当前 点 去 数据 集中 ， 种 子 点 己 可 随机 选 
取 ， 由 于 聚 类 针对 的 是 去 除 主干 后 的 点 云 数据 ， 
各 个 叶片 已 有 一 定 距 离间 隔 ， 因 此 ， 只 需要 选择 
合适 的 距离 阐 值 和 迭代 次 数 ， 聚 类 后 的 叶片 点 集 
数量 均 保 持 不 变 。 


3.2 棉花 叶片 对 应 关系 建立 


利用 2.3 节 介绍 的 方法 分 别 实现 不 同 生育 期 
作物 点 云 数据 及 叶片 间 的 一 一 对 应 关系 。 图 4 是 
1 号 棉花 第 46 天 和 第 49 天 叶片 中 心 三 维 欧 式 距 离 
构建 的 和 矩阵。 经 过 使 用 匈牙利 算法 计算 后 得 到 最 
佳 匹配 ， 阴 影 框 内 数值 表示 最 佳 匹配 时 对 应 叶片 
中 心 的 三 维 欧式 距离 ， 棉 花 植 株 7,, ,时 期 9 号 叶 
片 在 7T, 时 期 无 叶片 与 其 对 应 ， 为 新 长 出 的 叶片 。 

多 时 序 的 棉花 植株 间 各 叶片 的 对 应 关系 建立 
结果 如 图 5 所 示 ， 不 同时 序 点 云 中 的 同一 叶片 采 
用 同一 颜色 演 染 。 可 以 注意 到 ， 顶 部 的 新 生 叶 片 
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图 4 1 号 棉花 第 46 天 和 第 49 天 叶片 中 心 点 计算 的 最 佳 匹配 


Fig. 4 The best match between leaves of No. 1 Cotton Plant on 


day 46" and day 49" 


未 能 在 前 期 数据 中 找到 对 应 叶片 ， 而 植株 底部 的 
红色 叶片 随 着 植株 的 长 大 成 为 老 叶 脱落 。 


3.9 表 型 参数 计算 及 动态 变 


3.3.1 株 高 和 体积 

通过 在 提取 主干 过 程 中 得 到 的 作物 主干 部 分 
来 测量 棉花 作物 的 株 高 。 分 别 获 取 主 干 的 最 低 点 
和 最 高 点 的 三 维 坐标 ， 两 者 之 间 的 欧 氏 距离 即 为 
植株 高 度 。 将 通过 点 云 计算 的 作物 表 型 参数 与 人 
工 使 用 卷 尺 测量 的 表 型 参数 进行 对 比 ， 从 而 评价 
通过 作物 点 云 数据 提取 作物 表 型 参数 的 准确 性 。 
精度 评估 通过 绝对 误差 、 相 对 误差 及 线性 回归 决 
定 系 数 R (Coefficient of Determination) 来 实现 。 
其 中 绝对 误差 可 描述 为 作物 点 云 数据 计算 的 表 型 
参数 值 与 人 工 测 量 参 数值 的 差 值 ， 相 对 误差 可 描 
述 为 作物 表 型 参数 的 绝对 误差 与 人 工 测 量 的 参数 
值 之 比 乘 以 100%， 如 公式 (2) 所 示 。 

E, = |Y - L/L x 100% (2) 

其 中 ，E 为 相对 误差 ; 了 为 点 云 数据 计算 的 
表 型 参数 值 ; LAAT AM EBACE . 

R? 用 于 度量 作物 点 云 数据 计算 的 表 型 参数 与 
人 工 测量 参数 的 相关 性 ， 公 式 (3) 所 示 。 


Sry Yh 
_ i=l i=l 


R (3) 


3 ery 
其 中 ,7 为 作物 植株 叶片 数 ，7 为 人 工 测量 
的 表 型 参数 值 ; 了 为 人 工 测量 表 型 参数 值 的 平均 
fü; 这 为 点 云 计算 表 型 参数 值 。 
包含 给 定 三 维 空间 中 点 集 的 最 小 凸 多 面体 被 
定义 为 这 些 点 的 凸 包 ， 凸 包 的 生成 和 三 维 空间 中 
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(a)1 号 棉花 
(a)n € (b) 叶 长 估 测 
cA T ^. A 
= eae (c) 叶 片 寺 装 (A) AK 
= 图 6 ATE AR et en AR te 
de Fig. 6 Estimation for cotton plant volume, leaf length and 
e leaf area 
ei F 一 图 7 为 对 株 高 、 叶 长 、 叶 宽 等 参数 进行 精度 
e 评定 ， 通 过 对 12 株 棉花 所 有 计算 出 的 株 高 (如 
号 棉花 = yie "p 
i eed BUR). MPA A RENDERE GIG 
注 :从 左 至 右 分别 为 叶片 对 应 后 的 作物 点 云 在 第 46 天、49 X 58 " 
s i 则 
= 天 .63 天 的 结果 图 示 。 同 一 叶片 用 同一 种 颜色 泻 染 he 
> 图 5 3 株 棉花 4 个 时 间 点 各 叶片 间 的 对 应 关系 了 高 度 的 相关 性 ， 株 高 、 叶 长 和 叶 宽 的 决定 系数 
x Fig. 5 Leaf correspondence of three cotton plants at four dif- 4) 8| 1.0000, 0.9964, 0.9993, Hy tit F 1.0. 
a ferent growth stages 该 结果 表明 通过 三 维 激光 扫描 技术 可 获取 作物 
c 的 高 精度 表 型 参数 。 主 要 由 于 : 1) 所 使 用 的 近 
E Shy A te UE RE S = E aM . . i 
= 点 云 的 分 布 紧 密 相 关 。 因 此 ， 最 小 包围 凸 包 党 用 景 激光 扫描 系统 扫描 精度 较 高 ， 可 达 0.03 mm; 


来 表示 三 维 空 间 中 作物 植株 的 体积 。 棉 花 体积 计 2) 与 基于 影像 的 二 维 表 型 量 测 方法 相 比 ， 三 维 
算 过 程 : 输入 原始 植株 点 云 数 据 ， 通 过 MAT- 激光 点 云 数据 多 了 一 个 维度 的 信息 ， 因 而 在 表 型 
LAB 获 得 棉花 植株 的 凸 包 ， 即 可 自动 计算 喇 包 的 。 参数 测量 上 更 具 优势 。 

体积 ， 如 图 6 (a) 所 示 。 


s00 140 180 0.9786x + 0.4236 
7 © : y= x r 
g 550 | y=09985x-1.6975 P 了 人 dd P g 160 R? = 0.9993 e° 
E wn P É | ET P 
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= .9 ux 80 
IK 400 "i | K 60 e 
" AE 
350 A Ua e? 
300 LE 20 9 
300 400 500 600 20 40 60 80 100 120 140 20 70 120 170 
人 工 测 量 值 /mm 人 工 测量 值 /mm 人 工 测量 值 /mm 
(a) 株 高 测量 结果 (b) 叶 长 测量 结果 (c) 叶 宽 测 量 结果 


图 7 株 高 \ 叶 长 和 叶 宽 测量 结果 


Fig. 7 Measurment results of plant height, leaf length and leaf width 
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3.3.2 ”时 长 . 叶 宽 和 叶 面 积 估算 

棉花 叶片 存在 一 定 弧 度 的 弯曲 。 因 此 ， 直 接 
利用 叶片 两 端点 的 欧式 距离 值 来 表 描 述 叶 长 并 不 
适用 。 本 研究 采用 一 种 逐渐 和 逼 近 的 方法 来 测量 叶 
KRKE, (BI rn PA A (X, Yy, 
Z) MB (X,, Ys, Zs), AB 两 端点 间 的 距离 为 
Ls， 取 AB 中 点 C (图 6 (b) )， 分 别 计算 点 A、 
C 和 点 B、C 之 间 的 距离 和 Lsco S Llet Lc, 
WHE Lis M Las HKE mE K E d, EEE 7s 和 


Li EL, - L| «d. 令 : 


L,-7L,*d (4) 

WNIK; 否则 分 别 对 AC 和 BC 继续 处 
理 。 用 B 替 换 AC 段 中 端点 标号 C，AC 上 段 更 新 为 
AB 上 段 。 该 过 程 一 直 壕 代 执 行 ， 直 到 距离 差 值 满 
足 阔 值 条 件 为 止 。 叶 宽 也 采用 逐渐 禹 近 法 完成 
计算 。 

对 棉花 植株 点 云 数据 进行 叶片 器 官 分 割 得 到 
每 个 独立 的 叶片 点 云 数据 ， 使 用 Geomasgic Stu- 
dio 2013 软件 对 叶片 点 云 数据 进行 封装 ， 获 得 叶 
片 的 网 格 数据 ， 之 后 计算 各 网 格 的 面积 之 和 ， 即 
为 叶片 的 面积 。 在 网 格 化 过 程 中 ， 需 要 对 网 格 数 
量 和 距离 值 进 行 约束 ， 如 果 网 格 数量 越 多 ， 则 目 
标 表面 越 精 细 ， 计 算出 的 面积 值 与 真实 值 越 接 
近 。 本 研究 设置 的 网 格 数量 为 2,500,000 个 ， 距 
离 值 为 0， 因 此 ， 能 够 保证 得 到 足够 多 的 网 格 数 
量 ， 也 即 较 高 精度 的 面积 估 值 。 封 装 后 的 叶片 如 
图 6 (c) 和 6 (d) 所 示 。 

3.3.3” 表 型 参数 动态 变化 估计 

棉花 植株 点 云 数据 获取 过 程 中 进行 了 多 次 的 
作物 点 云 数据 采集 ， 因 此 可 以 将 同一 作物 不 同 生 
育 期 提取 的 表 型 参数 按照 先后 顺序 进行 对 比分 
析 ， 从 而 实现 作物 相 邻 时 期 表 型 参数 变化 过 程 ， 
并 将 其 结果 可 视 化 ， 为 育种 和 遗传 学 相关 研究 提 
供 数据 基础 。 

3 株 棉花 的 4 个 时 间 序 列 的 株 高 、 体 积 的 变 
化 趋势 如 图 8 (a) 和 8 (b) 所 示 ， 叶 长 、 叶 宽 、 
叶 面 积 的 动态 变换 趋势 分 别 如 图 8 (c) ~8 (k) 


化 趋势 较为 平缓 ， 呈 稳步 上 升 态势 。 体 积 的 增长 
趋势 和 株 高 较 相 似 ， 主 要 原因 在 于 体积 的 计算 和 
株 高 相关 性 较 强 。 老 叶 面 积 变化 不 大 ， 新 叶 面 积 
增加 趋势 明显 ， 如 1 号 棉花 中 顶部 的 三 片 叶子 
(图 5 (a) ) Porat Ben Re. REE., EK 
色 。 叶 面积 和 叶 长 的 变化 趋势 较为 一 致 ， 同 样 因 
为 叶 面 积 与 叶 长 有 较 强 的 相关 性 。 

本 人 研究 中 ， 从 1 株 棉 花 植株 的 点 云 数 据 采 
集 开 始 到 输出 各 表 型 参数 结果 值 ， 累 计 耗 时 约 
30 min， 其 中 数据 采集 约 15~20 min, EF HER 
和 叶片 聚 类 在 30 s 内 完成 ， 体 积 计算 约 11s， 单 
叶 的 叶 面 积 估算 约 耗 时 30 s。 此 过 程 中 ， 数 据 采 
集 较为 耗 时 ， 该 过 程 与 采集 对 象 的 复杂 度 有 关 ， 
对 象 越 复杂 ， 则 耗 时 越 长 。 


4 结论 与 讨论 


为 挖掘 三 维 激 光 扫 描 技术 在 作物 表 型 组 学 中 
的 应 用 ， 本 研究 利用 激光 扫描 技术 获取 了 棉花 的 
多 时 序 三 维 激光 点 云 数 据 ， 实 现 了 三 维 空间 内 株 
高 、 叶 长 、 叶 宽 、 叶 面积 及 体积 等 表 型 参数 获取 
和 动态 量化 。 主 要 取得 的 结果 如 下 。 

(1) 利用 RANSAC 算法 结合 直线 模型 提取 
了 3 株 棉花 共 12 组 点 云 的 主干 提取 ， 并 对 剩余 点 
云 使 用 欧式 聚 类 算法 完成 每 一 株 棉 花 的 叶片 分 
割 。 完 成 了 棉花 植株 的 叶 长 、 叶 宽 、 叶 面积 、 株 
高 和 体积 等 参数 的 计算 ， 对 比试 验 结果 表明 ， 提 
取 的 株 高 、 叶 长 、 叶 宽 与 人 工 测 量 值 均 具 有 高 度 
的 相关 性 ， 决 定 系数 R: 趋 近 于 1.0， 克 服 了 传统 
测量 方法 主观 性 强 、 具 有 一 定 的 破坏 性 等 缺点 ， 
为 后 续 产量 等 研究 葛 定 基础 。 

(2) 获取 同一 株 棉 花 作 物 不 同 生育 期 的 多 时 
序 点 云 数据 ， 研 究 了 作物 点 云 的 粗 配 准 和 叶片 器 
官 的 一 一 对 应 方法 ， 实 现 了 对 应 关系 的 可 视 化 。 
并 根据 吉 官 间 的 对 应 关系 实现 了 棉花 表 型 参数 变 
化 的 动态 量化 过 程 。 通 过 研究 茎 秆 、 叶 片 不 同时 
期 的 动态 变化 ， 可 为 剖析 棉花 个 体 的 器 官 生长 发 
育 过 程 提供 精确 的 信息 ， 也 可 为 构建 作物 生长 的 


所 示 。 从 图 中 可 以 看 出 ， 该 观测 期 内 ， 株 高 的 变 


动态 模型 提供 可 靠 的 数据 基础 。 
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图 8 3 株 棉 花 各 表 型 参数 随时 间 变 化 趋势 图 


Fig. 8 Dynamic change of phenotypic traits of the three cotons 
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Cotton Phenotypic Trait Extraction Using Multi-Temporal 
Laser Point Clouds 
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ty, Wuhan 430070, China; 4. College of Plant Science & Technology, Huazhong Agricultural University, Wuhan 430070, China) 


Abstract: To cope with the challenges posed by the rapid growth of world population and global environmental changes, schol- 
ars should employ genetic and phenotypic analyses to breed crop varieties with improved responses to limited resource environ- 
ments and soil conditions to increase crop yield and quality. Therefore, the efficient, accurate, and non-destructive measurement 
of crop phenotypic traits, and the dynamic quantification of phenotypic traits are urgently needed for crop phenotypic research, 
and breeding as well as for modern agricultural development. In this study, cotton plants were taken as research objects, and the 
multi-temporal point cloud data of cotton plants were collected by using three-dimensional laser scanning technology. The multi- 
temporal point clouds of three cotton plants at four time points were collected. First, RANSAC algorithm was implemented for 
main stem extraction on the original point cloud data of cotton plants, then region growing based clustering was carried out for 
leaf segmentation. Plant height was estimated by calculating the end points of the segmented main stem. Leaf length and width 
measurements were conducted on the segmented leaf parts. In addition, the volume was also estimated through the convex hull 
of the original point cloud of plant cotton. Then, multi-temporal point clouds of plants were registered, and organ correspon- 
dence was constructed with the Hungarian method. Finally, dynamic quantification of phenotypic traits including plant volume, 
plant height, leaf length, leaf width, and leaf area were calculated and analyzed. The overall performance of the approaches 
achieved a matching rate through a series of experiments, and the traits extracted by using of point cloud showed high correla- 
tion with the manually measured ones. The relative error between plant height and manual measurement results did not exceed 
1.0%. The estimated leaf length and width on point clouds were highly correlated with the manually measured ones, and the co- 
efficient of determination was nearly 1.0. The proposed 3D phenotyping methodology can be introduced and used to other crops 
for phenotyping. 


Key words: cotton phenotypic traits; LIDAR; stem extraction; leaf segmentation; point cloud registration; 3D phenotyping 
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